Sonoluminescing Gas Bubbles 



I. Scott 1 , H.-Th. Elze 12 , T. Kodama 2 and J. Rafelski 1 

1 Department of Physics, University of Arizona, Tucson, AZ 85721 

and 

2 Universidade Federal do Rio de Janeiro, Instituto de Ffsica 
Caixa Postal 68.528, 21945-970 Rio de Janeiro, RJ, Brazil 

Revised August 1998 

Abstract 

We draw attention to the fact that the popular but unproven hypothesis of shock-driven 
sonoluminescence is incompatible with the reported synchronicity of the single bubble sono- 
luminescence (SBSL) phenomenon. Moreover, it is not a necessary requirement, since we 
show that the sub-shock dynamic heating in gas bubble cavitation can lead to conditions 
required to generate intense lOOps light pulses. To wit we study the dynamics of the in- 
terior of a cavitating gas bubble subject to conditions suitable for sonoluminescence. We 
explore quantitatively the transfer of energy from the sound wave to the bubble interior, the 
frequency of atomic collisions in the bubble, the limits of quasi-stability of the non-linear 
bubble oscillations driven by an acoustical field, and obtain the implied reaction time scales. 

PACS: 47.40.Dc, 47.10+g, 43.30.+m, 78.60.Mq 
1. Introduction 

Current interest in the gas bubble cavitation dynamics has been revived by the development of 
experimental conditions which allow a single bubble 0, trapped at the anti-node of a stand- 
ing acoustic pressure wave, to emit relatively intense short light pulses 0. This single bubble 
sonoluminescence (SBSL) displays considerable sensitivity to the experimental conditions and pa- 
rameters. For a summary of the many ensuing experimental results we refer to the recent review 
by Barber et. al. ||. In particular they have found that the light flashes occur within 0.5 ns of 
maximum compression, last often less than 50 ps, comprise up to a million photons, and occur 
with synchronicity exceeding the frequency stability of the driving acoustical wave. The interplay 
of the diffusion processes and dynamical motion of the bubble has been also explored in depth [|j] 
suggesting that there is ongoing diffusion of the gas between the bubble and the liquid, and that 
the gas contents of the SBSL bubble is result of self-fine-tuning dynamics (rectified diffusion). 

The dynamics of the interior of the gas bubble and the energy inflow from the surrounding fluid 
into the bubble need to be understood in order to describe the energy focusing process in SBSL 
which is capable to convert acoustical energy into light. For sub-shock, spherical bubble dynamics 
we develop here a modification of the original Rayleigh-Plesset (RP) cavitation model ||, which 
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addressed only the dynamics of the spherical gas-liquid interface. The bubble interior enters the 
RP-model in terms of the (average) bubble gas pressure contributing to the surface dynamics and 
relevant in particular in the collapse phase. Our extension, based on variational approximation 
of hydrodynamics |6|, IT], accounts also for the kinetic energy of the gas in the bubble. In this 
approach we develop further the observation made earlier regarding the homologous bubble shape 
||, but we allow for an adiabatic change of the bubble density function between the collapse and 
expansion phases. Our sole present objective is to establish that there is rapid transfer of energy in 
sufficient quantity to the bubble, and that the bubble energy quenching occurs equally rapidly on 
sub-nanosecond scale. In this way we show that the required energy focusing can be accomplished 
just in a short time interval without inception of shock waves in the bubble, and that it arises from 
a combination of highly non-linear bubble dynamics and rapid quenching by acoustical radiation 
of bubble energy content |§ . 

Following a general discussion of the shock wave mechanism in the SBSL phenomenon in section 
2, we develop in section 3 the variational hydrodynamic description of the gas in the interior of 
the bubble. In section 4 we present numerical solutions and we study the dynamics of the energy 
flow to the gas bubble, and explore the dynamics of matter and energy within the bubble. We 
summarize our results in section 5. 



2. Shock waves and light pulses in SBSL 

Perhaps the most popular proposed SBSL explanation up to this day involves the hypothesis that in 
order to generate the ultra-short 50ps light pulses seen in SBSL |3|, [UJ one has to resort to extreme 
conditions that require shock wave formation in the bubble. Numerical hydrodynamic flow models 



incorporating inward shock formation [11, 12], as well as possibly outward motion shocks [13] have 



been studied. Despite considerable effort it is hard to prove convincingly that shock waves are the 



power-source behind SBSL ||14|| . This is so, since there is no clear understanding how the shock 
motion energy is converted into light pulses that depend very sensitively on the atomic composition 
of the bubble. The only reason we see for shock mechanism popularity is that the extreme matter 
conditions associated here with SBSL naturally last for a relatively short time period, well within 
the lifespan of the SBSL light pulse, which were originally reported to be less than 50ps long ||. 
More recently, depending on the details of the SBSL conditions, pulse-length in excess of lOOps 
has been reported [[K] . 

Since we can obtain such time scale in section 4 without shocks, one must ask if the shock 
hypothesis should not be abandoned, for there is a serious internal inconsistency: only perfectly 
spherical inward shock could be reflected without destroying the bubble. However, light emission 



asymmetry has now been observed [15|: there is a class of SBSL events in which significant light 
ellipticity in the bubble is required, at the level of 20%. It is natural to associate this asymmetry 
with matter distribution asymmetry. But then an inward shock wave could not bounce, the bubble 
matter would spray outwards, the bubble would burst. Even if a bubble were to be rapidly (i.e. 
within a cycle of the driving harmonic pressure field) recreated, there remains a contradiction with 
the well advertised magic of SBSL, that synchronous light pulses are emitted over billions of shock 
cycles with timing precision exceeding significantly that of the driving acoustical field generator. 
We recall that for this reason SBSL has even been proposed as a possible cheap and precise clock 
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Should we accept as argued that inbound shocks destroy the bubble, they will also erase bubble 
'memory'. The enhancement of temporal precision of cycle definition would have to arise from the 
liquid alone, a highly improbable situation, for which a mechanism is not in sight. An outward 
shock would not need to be reflected, but then the light pulse would be produced at the interface 
between gas and liquid, which is also inconsistent with the experimental results |L5| . In our opinion 
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the precision of light pulse periodicity can only arise from nonlinear dynamics within the gas bubble 
to which the external pressure field in the liquid is merely an energy source and hence it can be 
oscillating at a less well determined frequency. 

But what mechanism can produce light, if there is no shock? It is natural to expect that 
the light emitted is product of dynamical processes in the bubble. Any microscopic light pro- 
duction mechanism which does not generate strong spectral lines is suitable, including radiation 
from neutral atom collisions with continuum emissions [[17]]. This collision induced emission (CIE) 



mechanism for generation of SBSL-light pulses appears most consistent with the reaction frame- 
work we are proposing here. CIE process relies on relatively frequent and sufficiently energetic 
atomic collisions in high density matter. To address qualitatively the viability of this mechanism, 
we establish below in more detail the frequency of two-body atomic collisions at relative atomic 
energy that would allow the light pulse formation. While we do not pursue further more involved 
atomic physics aspects of the CIE process, we note that the densities reached in the bubble at 
sonoluminescence implicate three and more body collisions. 

We note that a bubble with radius similar to light wavelength is not opaque to radiation, 
which thus is not emitted from a surface. Light emitted from the bubble center will not show 
features characteristic of being in equilibrium with the (thermal) source. Therefore the conclusion 
that thermal radiation is excluded fl8f , is not at all providing evidence for shock- waves as the 



mechanism of SBSL. 

Another non-shock process has been proposed as mechanism for SBSL, involving strong, jet-like 
deformations invoked by the translational motion of the bubble in the gradient of the acoustical 



field [jnj. While jets indeed are likely to develop as the bubble begins to collapse, in our opinion 
the long term bubble stability implies here that the increasing internal bubble density and the 
increasing surface tension are capable to stabilize the spherical shape of the bubble preceding the 
short final collapse phase, during which the SBSL light is emitted f2~0fl . Were this not the case, the 



jet would have to excite SBSL in the final stages of bubble collapse with a non-negligible chance 
for bubble destruction. 

3. Variational bubble hydrodynamics 

Our theoretical extension of the standard RP approach concerns the dynamics of an 'adiabatically' 
homologous gas flow. This means that the shape and velocity of the interior matter are governed 
by the rapidly evolving radius parameter R(t) and surface velocity R(t), allowing for slower change 
in the homologous shape of matter characterized by a collective shape parameter a(t). In order 
to exploit this to the fullest, we implement the dynamics within the variational hydrodynamical 
method |7J], which has been discovered nearly 40 years ago (see Ref. ||, section 13, and references 
therein). Using the variational approach we sidestep the need for a full hydrodynamic solution 
of the gas motion. This simplified approach is suitable only for exactly spherical, (quasi) one- 
dimensional motion, and corresponds to a variational solution of the hydrodynamic problem with 
irrotational flow. 

We consider a gas bubble of radius R(t), shape g(x,a(t)) mass Mq = ^nR^po (expressed in 
terms of the initial radius Ro and density p ) surrounded by an incompressible liquid of density 
Pl- in terms of the -R-scaled radial variable, the mass-conservation constraint, Mq = 0, implies: 

M G r 

pG = ^RW 9Mt)) ^ X = W) <1, (1) 

[ 1 x 2 gdx = l, [ 1 x 2 ^-dx = 0. (2) 

Jo Jo oa 

R(t) and the collective variable a(t) are the dynamical functions we will consider. a(t) comprises 
two important physical phenomena: 
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1. the slow underlying variation ('drift') describes how the shape of the mass density in the 
bubble changes between the contraction and expansion phases of the bubble motion - solving 
at fixed a = Const. We are implying that the shape shift is adiabatic compared to the rapid 
matter flow, expressed by changes in R; 

2. small, but rapid oscillation of a(t) may be superposed with this slow drift; they correspond 
to the acoustical cross-talk oscillations between the surface and the interior of the bubble 
with frequency uj ~ c s /R < 10 9 s _1 , where c s is the sound velocity in the gas. 

The introduction of a(t) rather than time t in definition of the shape function g in Eq. ([!]) serves 
to separate these two types of non-homologous dynamics of the bubble. 

In order to determine the dynamics of the gas bubble, we will seek an extremum of the action 
for the combined gas-liquid system using as dynamical variables R(t), R(t) a(t), a(t) . The matter 
shape function g(x, a) is imposed by the external constraints. The kinetic energy of the gas can 
be easily expressed as a function of g: solving the continuity equation we obtain the local flow 
velocity in the gas bubble: 



X 



v(r) =xR-Rad } 5 = / dx'x' 2 ^- — > (3) 

x l g[x\ a) Jo oa 

After some integrations by parts of the time integral (action) the contribution in action is found 
= 

K G = ^{-I lR R + I 3 a 2 R 2 } , (4) 

with: 

h — dxx 4 g, h = dxx 2 g5 2 . (5) 
Jo Jo 

When 8 is neg lected (§2^0) one is in the exact homologous limit. 

Eq. (^]) complements the kinetic energy of the incompressible liquid of density p^. Solving the 
continuity equation we obtain the local flow velocity and thus also the kinetic energy in the liquid: 

v L = R^, r>R; K h = Aiip h R^. (6) 
The internal (potential) energy of the gas is: 

W G = W G [R,g(x;a)} = J e(r,t) d 3 r = M G £ dxx 2 gS , (7) 
where e is the energy density and 

£ = E/N = e/p 

is the specific energy of the gas. The total potential energy W of the bubble (normalized to zero 
for a bubble at rest under ambient conditions (pressure P = Pq and temperature T = T ) is: 

W = W G + 4vra( J R 2 - R 2 ) - P L (V - V G ) • (8) 

We note that the second term in Eq. (§) is the energy from the fluid-gas interface, and Pl = 
-Po + Pa(t) is the pressure in the liquid arising from the atmospheric pressure and any further 
time dependent externally applied (acoustical) pressure. Setting Ji = J 3 = = Kg i.e. ignoring 
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the interior dynamics of the bubble, and varying Kl — W with respect to R(t) one obtains the 
(improved) RP-model equations [|2T| . 

We first determine the dynamical equations which characterize the drift in a(t). The tool of 
convenience here is the Hamiltonian for the motion of this dynamical variable. The generalized 
momentum following from Eq. (U) is: 

n a = ^ = dI 3 (a)R 2 M G , (9) 
da 

and the Hamiltonian for the gas is then 



n G (a,u a ) = ^ a + v a . (10) 

Where the (collective) mass Ai a and potential V a are: 

M a = I 3 (a)R 2 M G , V a = ^h(a)RR + M G f dxx 2 g(a)S[p(a)} . (11) 

2 Jo 

where I\, I 3 are given in Eq. (|5|). The rapid oscillations inherent in the harmonic Hamiltonian 
presented are centered around the minimum of V a , which drifts as the values of R, R change. At 
the same time there is a drift in the mass Ai a which impacts the frequency of the rapid oscillations, 
but not the location of the minimum of V - thus if we study solely the drift in a-coordinate, we 
are effectively neglecting the deviations from the homologous shape inherent in the 5-dependence 
ofJ 3 . 

As we search for the minimum of the effective potential in order to derive equations governing 
this motion we will need to invoke a few simple thermodynamic properties. From the first law of 
thermodynamics we recall the relation of pressure P to the internal specific energy 8: 

p — m = P *m . (12) 

The heat function (enthalpy) h is introduced in view of the relation: 

d ,„ P, . ldP . „ P r dP 

— (£ + -) =-— — > h = £+-= —. (13) 

dp p p dp p JS=Const p 

The minimum of V a determines the drifting value of a , which is obtained by solving: 

1 f'^idg r rj u , 1 



M G da 



J dxx 2 -^\h[p\ -h c + -x 2 RR^ = 0. (14) 



Here we have introduced a constant h c , which integrates to zero by the second expression in Eq. @, 
and which assures that the solution of Eq. ( |T4] ) is found where the curly bracket in Eq. ( II ) vanishes. 
Considering the limit x 2 —>■ we see that h c = h(x = 0). Differentiating with respect to x 2 the 
argument of the curly bracket we obtain, under adiabatic constraints and using the definition of 
h, see Eq. (fTJ): 

-2 RR= dx^ {h[p] - K}= dx^WdP> (15) 



and thus: 



d . . RR 9 1 \ 9P . . 

^ 9{x ' a) = -2^MY c ^ (p) = ^' (16) 



6 



Scott, Elze, Kodama and Rafelski 



where c s is the local sound velocity in the gas. Here the factor RR is to be seen as a given 
parameter; in reality it is determined by the usual RP-dynamics, which does not depend much on 
the inclusion in the dynamics of the bubble interior kinetic energy of the. Thus Eq. ( pf) can be 
solved independently of the RP-dynamics, yielding the interior shape of the bubble at each phase 
of the dynamical evolution of bubble radius. We note that the sign of the RP-force factor RR in 
Eq. flnS| ) determines the gradient of the density shape function: when the bubble is collapsing, the 
surface density is higher than the interior density, while when it is expanding the opposite is the 
case, and a flat distribution is found in the force free instants of the bubble motion, when R = 0. 
We have explored elsewhere other theoretical features not discussed in depth here, such as the 
small fluctuations in a(t) f7j, the back-reaction of the bubble dynamics into the RP-equations @, 
non-perturbative treatment of the sound emission^. Also, we have obtained Eq. (fl6|) for g(x,a) 
directly from the Euler equation in the limit 5 — > . 

We note that only for a particular choice of the initial condition g(x=0) can we obtain a properly 
normalized solution of the nonlinear differential equation, Eq. (|T6|) . However, this we can do easily 
for any EOS subject to an adiabatic constraint, and specifically also for the vanderWaals system 
we consider here, with finite size correction to the volume in the polytropic equation of state: 



Pm ,V ( P V Lr , IP 



P[p] = Po ~ 1 



7 - — 

Pm 



e\p] 



p 



7 - 1 \p p 



PO j \Pm - PJ 7 - 1 P 

(17) 

with p m = M G /a 3 = SMq/Attt^ and r m is the vanderWaals excluded radius. 



4. Physical properties of sonoluminescing bubbles 

The motion of the bubble is primarily damped by sound radiation in the crash phase. We have 
determined that in rough detail a full treatment of the sound radiation is equivalent to first order 
in R/c description when adjustment (i.e. fit) of the damping constants is allowed for ||. We thus 



follow the approach of Lofstedt |2T| and account for the sound damping influence on dynamics 
of R(t) by including the appropriate first order in R/c terms in our extended Rayleigh-Plesset 
equation. We also incorporate the damping due to liquid viscosity. We do not study the effect 
of compressibility of the fluid surrounding the gas bubble. We use the parameters adjusted in 
Ref. [^1]] to fit the dynamics of the bubble radius. 

The standard example we consider here is an argon gas bubble comprising N G = 10 10 atoms 
immersed in water. The parameters of the liquid are: density pl — lg/cc, surface tension a = 
0.03 kg/s 2 , sound velocity ch 2 o — 1450 m/s, and viscosity rj = 3 • 10~ 3 kg/(ms) . The equilibrium 
radius of a bubble under ambient pressure P$ = 1 atm at ambient temperature T = 277K obtained 
from the vanderWaals equation state is Rq = 4.39pm and the excluded radius parameter is 
a = 0.503 pm . The acoustical driving pressure is taken to be P&(t) = P a sin(a;t), with P a = 1.4 atm 
and uj/2tt = 25 kHz . 

Our numerical procedure is straightforward: we determine the solutions of the RP-equations, 
and use these to solve Eq. (|1(J) at a given instant in time. The solutions turn out to be very 
smooth in time and space and we can interpolate quite easily the solutions with large time steps. 
Subsequently, we can study any property of the bubble. 

When we evaluate, as a function of time, the energy content in the bubble, we find that there 
are three time scales of interest in the gas bubble dynamics: the overall energy of the liquid and 
the bubble varies on scale of driving frequency, within the collapse window of a few microseconds, 
as can be seen in Fig. la. Note that we present the energy of the liquid per atom in the gas bubble. 
The solid line is the total energy, the dotted line is the kinetic energy of flow, the dashed line is the 
internal (potential) energy arising from the work done on the liquid. Both components contribute 
significantly to the overall energy of the liquid. We see that the potential energy dips during the 
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slow expansion phase of the bubble motion, and then smoothly rises to a maximum value at which 
point it declines at a slightly faster rate. During the collapse the kinetic energy of the fluid rises to 
O(60 eV) per atom in the bubble gas, which energy is rapidly lost to elastic and inelastic processes 
when the bubble hits the van der Waals hard core. The dynamical time scale of the fluid, defined as 
the time difference in minimum-maximum-minimum energy evolution, is (9(10 /is), corresponding 
to a quarter of the acoustical driving cycle of 40/is. The interior gas energy varies on the scale of 
nanosecond, as can be seen in Fig. lb. In this example about 7% of the energy is transferred from 
the liquid to the gas. Finally, in Fig. lc we show a time-resolved view of the energy contents of the 
interior of the bubble. We see that most of the energy rise occurs within a 100 picosecond time 
window. The energy of the bubble disappears rather slowly; there is no viscous damping included 
in the dynamics of the gas. 

We now turn our attention to estimate the intensity of the sonoluminescence arising from the 
atomic collision processes; for this we obtain the collision frequency of all atoms in the bubble: 

f N = Na^-v Tel . (18) 
The effective collision cross section is assumed to be geometric and thus a ~ (47r/3a 3 /iV) 2 / 3 . A 



measure of relative velocity is obtained from the total energy per atom t> re i — y2E/m. We can 
obtain these results for a number of different conditions of driving pressure, ambient pressure, 
frequency, etc. In Fig. 2 we present the results for P a = 1.4, 1.3, 1.2 atm, selecting an arbitrary 
time axis origin so that these results can be shown together. The number of collisions shown is 
in 10 10 /ps, thus during the light flash duration of 50 ps we expect about 10 12 -10 13 atom-atom 
collisions to occur in the bubble. We recall that experimentally observed is the emission of up to 
million photons in the spectral range of visible light during the same period of time (see e.g. legend 
of figure 1 in Ref.0.) We also see in Fig. 2 that the number of collisions is changing smoothly, 
and in a manner corresponding to the emission intensity observed in the experiment [Q. The 
vertical lines indicate when in time the average energy per atom exceeds 1 eV. The realm of the 
frequent collisions at 1 eV and above is well within the 500ps time interval of the most extreme 
density, as is experimentally reported for SBSL ||. In our opinion these results strongly support 
the hypothesis that atomic collisions are the source of the short SBSL light flashes [ 17| , since the 



number of atomic collisions capable to generate visible photons exceed by the comfortable margin 
of 10 6 -10 7 the photon yield. 

While the onset of SBSL is a gradual process, depending on the magnitude of the driving 
acoustical pressure field, and could be explained as the onset of frequent high energy atomic 
collisions, the disappearance of SBSL when P a > 1.4 atm is quite sudden. We considered the local 
sound velocity c s (x) in the bubble and compare it to the local flow velocity v, Eq. (§), and have 
determined that this condition corresponds to the interior flow velocity in the gas bubble slightly 
exceeding, near to the surface, the local sound velocity. The flow velocity in the bubble, in units 
of local sound velocity, is shown for several values of driving pressure P a = 1.4, 1.3, 1.2 atm in 
Fig. 3 As P a — > 1.4 atm, we see in Fig. 3 that the local flows in the bubble are bound to occur with 
supersonic velocity. For the development of shock waves just this supersonic condition is required. 
Since in the experiment P a ~ 1.4 atm is the upper limit for the SBSL process, we are led to the 
hypothesis that it is the onset of shock instability in the bubble, which is responsible for the upper 
limit in P a for which SBSL can occur. 



5. Discussion and conclusions 



In summary, we have developed a semi-analytical dynamical model to describe the properties of 
the interior of the cavitating gas bubble. In our approach we have simplified the complexity of 
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the hydrodynamic gas flow problem by exploiting in a variational approach the nearly homolo- 
gous bubble dynamics. Considering that SBSL occurs in an extremely small domain of a large 
parameter space, our approach is in some aspects more capable to gain physical insights, than is 
full hydrodynamics. Several approximations were made: we did not address quantitatively in our 
numerical work the sound fluctuations inside the bubble during its motion. We did not consider 
here the back-reaction from the dynamics of the gas bubble to the dynamics of the liquid which, 
while important in principle, does not change the general dynamic behavior of the bubble . Nei- 
ther did we include the compressibility in the description of the surrounding fluid motion, but our 
state of the art implementation of the dynamics of the bubble includes in a way proposed in Ref. 

the dissipative effects such as viscosity of the fluid, and most importantly, acoustic radiation 
damping 21]. We checked and agree with Ref. |nj that allowing for compressibility of the 



liquid in the extreme conditions of SBSL further sharpens the characteristic time scale involved, 
as the denser liquid can more rapidly acoustically dissipate the air bubble energy. We find that 
even without this effect there is a natural 100-picosecond time scale in argon cavitation in water, 
arising from the transfer of kinetic flow energy from the liquid to the gas bubble in its expansion 
phase, followed by inertial collapse of the energized bubble. 

We were able to derive the density shapes of the bubble and to evaluate the energy content 
of the bubble as function of time. We have shown that the bubble energy has a characteristic 
relaxation time (9(1 ns), while the scale (9(100 ps) describes the rise in energy of the bubble. We 
have shown in quantitative fashion that it is likely that atomic collisions are the source of the light 
flashes seen in SBSL. We refrain here from the evaluation of the CIE rates [17]], pending further 



clarification of microscopic mechanisms (specifically the role of noble gases, three body collisions), 
and also since this task transcends the scope of this research project. Thus while this is not the 
time and place to explore in detail the light spectrum that can be produced, we note that during 
the short time interval allowing deep atomic collisions, CIE will produce light pulses of similar 
temporal length for different wavelengths in the visible range. 

The dynamical properties of the RP-model lead us to believe that it is the (inward) shock insta- 
bility that poses an upper-applied-acoustical-pressure-limit on the SBSL phenomenon. Specifically, 
when the strength of the driving acoustical pressure reaches the value at which experimentally 
SBSL phenomenon suddenly disappears 0, our calculations show that the internal bubble dy- 
namics must evolve faster than local sound velocity. This should lead to flow instabilities, which in 
all likelihood can destruct the bubble. Our results imply that the complex, near-shock-nonlinear 
dynamics of the bubble interior combined with light emissions in atomic collision processes can 
provide a suitable framework to explain the SBSL phenomenon. 
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Figure 1: Argon bubble with 10 10 atoms in water. Evolution of gas bubble energy: dotted line is 
the kinetic energy, dashed line is the potential energy, and solid line is the total energy. The thin 
solid line over arbitrary scale shows the sum of ambient pressure, 1 atm, and driving acoustical 
pressure (cu = 25 kHz), amplitude at bubble surface 1.4 atm: a) The energy in the fluid in eV per 
atom in the bubble; b) The energy in the bubble over nanosecond time scale; c) The fine resolution 
in time of the bubble energy during most extreme collapse. 
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Figure 2: Number of atom-atom collisions (times 10 10 ) per picosecond in the bubble as function of 
time for three values of driving acoustical amplitude, P a = 1.2, 1.3, 1.4 atm. Dashed vertical lines 
delimit range for 1 eV energy per atom in the bubble. 



12 



Scott, Elze, Kodama and Rafelski 



1.4 




Figure 3: Local flow velocity in units of the local sound velocity (termed Mach as function of 
the scaled radial coordinate of the bubble, r/R. Dashed line: Mach #=1.0. Results for driving 
pressure P a = 1.4, 1.3, 1.2 atm are shown (top to bottom). Event time chosen to maximize the 
local Mach # at the bubble surface. 



